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ABSTRACT 

Bright, ultracompact X-ray binaries observed in dense star clusters, such as Galactic globular clusters, must have formed 
relatively recently, since their lifetimes as persistent bright sources are short (e.g., ~ 10 8 yr above 10 36 ergs _1 for a 1.4M Q 
neutron star accreting from a degenerate helium companion with an initial mass of ~ O.2M ). Therefore, we can use the present 
conditions in a cluster core to study possible dynamical formation processes for these sources. Here we show that direct physical 
collisions between neutron stars and red giants can provide a sufficient formation rate to explain the observed numbers of bright 
sources. These collisions produce tight, eccentric neutron star - white dwarf binaries that decay to contact by gravitational 
radiation on timescales ~ 10 6 - 10 I0 yr, usually shorter and often much shorter than the cluster age. 

Subject headings: binaries: close — galaxies: star clusters — globular clusters: general — hydrodynamics — 
stellar dynamics — X-rays: binaries 



1. INTRODUCTION 

Ultracompact X-ray binaries (UCXBs) are persistent, bright 
X-ray sources (L x ~ 10 36 - 10 39 ergs _1 ) containing a neutron 
star (NS) accreting from a low-mass, degenerate compan- 
ion in a very tight orbit of period P < lhr. UCXBs may 
well be dominant among the bright low-mass X-ray bina- 
ries (LMX Bs) observed in old globular clusters (GCs), both 
Galactic (Deut sch et alJl2000t Ivan der Sluvs et al.H2005l) and 
extragalactic (Bildsten & Delove 2004). It was recognized 30 
years ago that the total numbers of LMXBs observed in GCs 
clearly indicate a dynamical origin, with formation rates ex- 
ceeding those in field populations by several orders of magni- 
tude (Clark 1975). Indeed, the stellar encounter rate in a clus- 
ter core is an excellent pr edictor for the presence of a bright 
LMXB (Jorda n et al 120041) . 

The growing importance of UCXBs is clear from the role 
they have played recently in a number of different contexts. 
They may dominate the bright e nd of the X-ray luminos - 
ity function in elliptical galaxies (Bildsten & Delove 2004). 
They pose a number of challenges to, and may allow us to 
test our fundamental physics of, stellar str ucture for low-mass 
degen erate or quasi-degenerate objects (iDelove & Bildstenl 

2003) . They may also connect in a fundamental way to 
NS recycling, as suggested by the fact that three out of 
six accretion-powered m illisecond X-ray pulsars known in 
our G alaxy are UCXBs (Chakrabarty 2004; Grebe nevet all 

2004) . Finally, UCXBs may well be the progenitors of the 
many eclipsing binary radio pulsar s with ver y low-mass com- 
panions observed in GCs (Rasio et al. 2000). 

Several possible formation processes for UCXBs are pos- 
sible. Exchange interactions between NSs and primordial 
binaries provide a natural wa y of f orming possible progen- 
itors of UCXBs llRasioetal.1120001) . This may well domi- 
nate the formation rate when integrated over the entire his- 
tory of a cluster. However, it is unlikely to be significant for 
bright UCXBs observed today. This is because the progen- 
itors must be intermediate-mass binaries, with the NS com- 
panion massive enough for the initial mass transfer (MT) to 
become dynamically unstable, leading to common-envelope 

1 Northwestern University, Dept of Physics and Astronomy, Evanston, IL 
60208. 

2 Dept of Physics and Astronomy, Vassar College, Poughkeepsie, NY 
12604. 



(CE) evolution and significant orbital decay. Instead, all 
main-sequence stars remaining today in a GC have masses 
low enough to lead to stable MT (and orbits that expand 
during MT, leading to LMXBs with wide periods and non- 
degenerate donors). Alternatively, some binaries with stable 
MT could evolve to ultra-short periods by magnetic braking 
dPvlvser & Savoniielll988HPodsiadlowski et al Jl2()'()2). How- 
ever, producing UCXBs through this type of evolution re- 
quires very careful tuning of initial conditions and is therefore 
v ery unlikely to e xplain most sources (?). 

Verbunt ( 1987) first proposed that a physical collision be- 
tween a NS and a red giant (RG) could lead to UCXB for- 
mation. In his scenario, the collision was assumed to lead 
directly to a CE system in which the NS and RG core would 
quickly inspiral. However, RG-NS collisions that occur now 
in old GCs (where RGs have low masses, close to m t0 ) do not 
lead to CE evolution. Instead, the RG envelope is promptly 
disrupted, leaving behind a n ecc entric NS-WD binary, as 
shown by Rasio & Shaping (119911) using 3-D hydrodynamic 
calculation. Nevertheless, if the post-collision NS-WD bi- 
naries can decay through gravitational-wave emi ssion all the 
way t o contact, they can still become UCXBs (iDavies et alJ 

2. OUTCOME OF COLLISIONS 

Using the 3-D Smoothed Particle Hydrodynamics (SPH) 
code StarCrash 3 we have computed about 40 represen- 
tative collisions betwe en various RG stars and a 1.4 M Q NS 
( Lombardi et alJ 120041) . In our models, both the NS and 
the RG core are represented b y point masses cou pled to the 
gas by (softened) gravity only (Rasio & Shapiro 1991). Our 
initial RG models were calcul ated using the stellar evolu - 
tion code described in de tail in [P odsiadlo wski et all (2002), 
Ivanovaetal. (2003), and Kalogera et al. (2004). The mod- 
els include stars on the subgiant branch (with total mass 
M = 0.8M Q , core mass m c =0.10M Q and radius R lg = l.6R Q 
and with M = 0.9M Q , m c = 0. 12M© and R Ig = 2R Q ), and sev- 
eral models near the base of the RG branch. Our most evolved 
models have M = 0.9M©, m c = 0.25M Q and R Ig = 6.8R Q . 
More evolved RGs contribute very little to the total colli- 
sion rate (see §3). The distance of closest approach for the 
initial collision varies from r p = 0.1 R lg (nearly head-on) to 

3 See http : / /www .astro. northwestern ■ edu/StarCrash/| . 
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FIG. 1. — Dependence of the gravitational radiation merger time on post- 
collision semimajor axis a and eccentricity e for a binary consisting of a 
1.4Mq NS and a 0.25 Mq WD. The points with different symbols show the 
results of our SPH calculations for six different giant models (starred sym- 
bols are for a 0.9 Mq star; others are for a 0.8 Mq star at different evolution- 
ary stages, e.g., triangles represent a subgiant, and circles our most evolved 
RG model! The symbol area is proportional to the collision rate, accord- 
ing to eq. J4j, i.e., symbols for less evolved RGs appear larger. The hatched 
area shows how the merger time (here for the line of constant merger time 
T gw = 10 lo yr) changes when we vary slightly the binary parameters: the up- 
per boundary corresponds to a 1 .5 Mq NS with a 0.45 Mq WD, and the lower 
boundary corresponds to a 1.3Mq NS with a 0.15Mq WD. 



r p = 1.3 R Ig (grazing). 

In agreement with previous SPH calculations 
jRasio & Shapiro! 119911 iDavies et all 119921) . we find that 
all collisions produce bound systems in which the RG core 
ends up in a high-eccentricity orbit around the NS. However, 
in contrast to those older studies, our new SPH calculations 
extend over much longer times (up to ~ 500 successive 
pericenter passages), allowing us to determine acc urately the 
final parameters of the orbit ( Lom bardi et all 2004). Typically 
~ 50% of the RG envelope is ejected to infinity, while most 
of the rest becomes bound to the NS. Only about ~ 0.1 Mq 
remains bound to the RG core, which will eventually cool to 
a degenerate WD (cf. §4). The material left bound to the NS 
will attempt to form an accretion disk as it cools. The fate 
of this material is rather uncertain. It could be accreted onto 
the NS and spin it up (in ~ 10 6 yr at the Eddington limit), 
or, more likely, it could be ejected if the energy released by 
accretion couples well to the gas. With an efficiency e, the 
entire mass of gas could be ejected to infinity in as little as 
T gas ~ 500(e/0.1) _1 yr. This very short lifetime justifies our 
assumption that the parameters of the post-collision orbits 
determined by our SPH calculations are nearly final, i.e., that 
the orbital parameters are no longer affected by coupling of 
the orbit to the residua l gas. 

When we apply the lPetersI Jl964) equations to these post- 
collision systems, we find that most of them inspiral on rather 
short timescales (Fig. 1). Therefore, we assume for the rest of 
this paper that all RG-NS collisions can produce UCXBs. 



3. COLLISION RATE 

Consider a NS of mass m ns in the core of a cluster contain- 
ing iV* ordinary stars (here we neglect binaries; see §4). If all 
these ordinary stars were turn-off stars of radius R to and mass 
m t0 , the collision rate for the NS would be 

K to = 2 7 rG(m to + m ns )A^ to cr- 1 V c - 1 , (1) 

where a is the relative velocity dispersion, V c is the core vol- 
ume. Here we assume that the collision cross section is dom- 
inated by gravitational focusing . 

To compute the collision rate with RGs, we take into ac- 
count that the number of RGs, dN„, within any small range 
of radii between R„ and R„ + dR Ig is proportional to the 
time dt spent there by the star as it ascends the RG branch, 
dN Ig = f Ig N*dt /t. Here f Ig is the fraction of stars with masses 
close enough to the turn-off mass to have become RGs, and r 
is the total lifetime (from the ZAMS to the end of the RG 
stage) of a turn-off star, only slightly larger than the clus- 
ter age. For a simple analytic estimate, we use the follow- 
ing approximate relation b etween age and radius (eq. A9 of 
iKalogera & WebbinU 19961) . 

/ fx" - 28 
R Ig (t)~R ZAM s[l--) , (2) 

where we use /?zams — 0.1R to - Next we replace dt by 
dR Ig /(dR Ig /dt) = (R to /R Ig ) 4 - 6 TdR Ig /R to . The collision rate 
for a NS with RGs between R Ig and R lg + dR lg is 

dTZ = 2.67rG(m t0 + m^R^a' 1 V C _I dN lg . (3) 

Here a collision is defined to be any encounter with a dis- 
tance of closest approach less than 1 .3 R Ig , consistent with our 
SPH results. Integrating this over R ls from the base of the RG 
branch, defined by setting R lg = bR to , to the maximum radius 
of a RG, R max ^ bRta, we find a total collision rate 

KvcxB^OSlf^b- 2 - 6 ^. (4) 

Alternatively, note that we could also directly integrate eq. (0 
over time, without changing the variable from t to R rg . Be- 
cause the collision rate is linearly proportional to radius when 
gravitational focusing dominates, we can then write 

7^-ucxb — 2.67rG(m to + OT ns )A^*/ 1 g^ rg cr" 1 V c " 1 , (5) 

where /? rg is the time-average radius of the RG. Using eq. 
it is easy to show that eqs. and (|5} agree. Eq. l|5} has the 
advantage that any stellar evolution treatment can be used to 
determine R lg , including fitting formulae more detailed than 
eq. or numerical results from stellar evolution calculations. 

The steep inverse dependence on b in eq. © indicates that 
the collision rate is completely dominated by the smallest 
RGs: although the cross section increases (linearly) with ra- 
dius, the faster stellar evolution at larger radii dominates, so 
that collisions are much more likely to happen when the star 
is just leavin g the main sequ ence, i.e., on or close to the sub- 
giant branch ( Verbun tll987l) . The corresponding core mass is 
also small, typically ;« c ~ 0.1 M Q for m to ~ 0.8 Mq. 

We now proceed to evaluate f rg . This depends on the 
mass function of stars in the cluster core, which we expect 
to be very different from the IMF because of mass segre- 
gation. Indeed, observations of cluster cores reveal flat o r 
even slightly rising mass functions (e.g., Richer et al. 2004). 
Here we assume that the number of stars within dm is pro- 
portional to m a dm, with a > -1, between a minimum m m \ n 
and a maximum m lo + Am. The spread Am of masses along 
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TABLE 1 

UCXB FORMATION IN GALACTIC CLUSTERS . 



Cluster 


frg/T 


/? rg 


logpo 


<7 


log Mi 


Nioo 


A'min 


NGC 1851 


0.071 


6 


5.7 


10.4 


6.0 


0.11 


0.85 


NGC 6624 


0.087 


4.2 


5.6 


5.4 


5.2 


0.14 


0.18 


NGC 6652 


0.076 


5.4 


4.8 


5.9 


5.4 


0.02 


0.05 


NGC 67 12 


0.070 


5.9 


3.0 


4.3 


5.0 


0.0005 


0.0004 


NGC 7078 


0.034 


7.1 


6.2 


12.0 


6.1 


0.16 


1.62 


Terzan 5 


0.10 


4.3 


6.1 


10.6 


5.6 


0.27 


0.87 


47Tuc 


0.081 


4.9 


5.1 


11.5 


6.1 


0.23 


0.23 



NOTE. — The RG lifetime fraction t K (r and the average RG radius R K (in Rq) 
are calculated directly from our stellar evolution code and used in eq. |5| (with a = 0) 
and eq. |5j. po is the cluster core density (in Mq pc~ 3 ), a is the ( 1 -D) velocity disper- 
sion (in kms~') and M t is the total cluster mass (in )■ Fo r Ter 5 and NGC 6652, 
logpo is based on the luminosity density from Diorsovski ' 1993) and an adopted 
mass-to-l ight ratio of 2. The value of log pp for N GC 6652 appears rather uncertain 
(see. e.g. JPrvor & MevlailT 993 : Diorgovski 1993). Values of logM for NGC 6652 
and Ter 5 are from Gnedin et al. 2002); a for Ter 5 is from Gnedin et al. 2002) 
and for NGC 6652 from Webbink ' 1985). Otherwise logpo, a and logM ( are from 
lPrvor&MevTanlll993ll . 



the RG branch is obtained from the mass-dependence of 
the main-sequence lifeti me f ms . Adopting the simple scal- 
ing f ms = T(m to /m) 3 ' 6 JHurlev et alJ |$QO0) we get Af ms = 
3.6r(m to /m) 3,6 Am/m. Setting Af ms = f rg , the total time spent 
on the RG branch, and m ~ m to gives Am ~ 0.28 m lo t rg / r. 

To be consistent with our previous definition of a RG hav- 
ing a radius > bR to and using again eq. (|2j we derive t rg /r = 
0.3fr~ 36 . We can now calculate / rg directly from IMF. Assum- 
ing Am < m to , and m min < m to we get 

f = ( a + 1)— ~ 0.28(a+ 1)^ = 0.08(a+ 1)£" 3 ' 6 . (6) 
m to t 

Combining this with eq (|4) we obtain the result, 

ft UC xB^0.04^ft to . (7) 

In steady state (justified given the short lifetimes ?ucxb *C t 
of the bright UCXB phase) we can then estimate the number 
of UCXBs per 100 NSs at present in a cluster as 

/Vioo — 1007?-ucxBfucxB- (8) 

The lifetime fucxB depends on the minimum luminosity for 
a system to be classified as an UCXB. For our estimates we 
adopt a minimum luminosity comparable with the observed 
minimum in our Galaxy, L x ~ 1 36 ergs~ 1 . The co rresponding 
lifetime is f UC xB ^ 10 8 y r (e.g- lRasio et all20 00). 

The present mass of a cluster M t is always less than its ini- 
tial mass Mt;o = f^l L M t , where /ml is the total mass loss frac- 
tion. About 40% of the initial mass is lost just through stellar 
winds and SN explosions, so that /ml < 0.6. Without tidal 
mass loss (Joshi et al. 2001 ), and adopting a lower mass cut- 
off of 0.1M© in the IMF of Kroupal (120021) . w e expect about 
1 NS p er 65 M Q of mass at present (see also llvanova et all 
l2004al) . About 5% of these NSs will be retained, depending 
on the escape velocity and th e NS natal kick velocity distri- 
bution (llvanova et alJl2004ah . The corresponding minimum 
number of UCXBs expected (without any tidal mass loss) is 
then 

7V min ~ 8 x 10- 4 M t ^ UC xB fucxB , (9) 

where M t is in Mq . 

In Table 1, we show numerical results for several Galac- 
tic clusters: all clusters where a UCXB has been identi- 
fied, and 47 Tuc (which does not contain any bright LMXB). 



The probability of finding a bright UCXB in a cluster like 
47 Tuc is only about 23%. NGC 6652 has poorly mea- 
sured parameters (see Note for Table 1), and our numbers 
for this cluster are necessarily uncertain. Two clusters, 
NGC 6624 and NGC 6712, are thought to have very eccen- 
tric orbits and to be on the verge of complete d isruption in the 
Galactic tidal field (iRichtler et al.lfT994l iGnedin etal]|1999t 
lAndreuzzi et al.l2001l) . This suggests that they may have had 
much higher mass and density in the past. Indeed, observa- 
tions show that NGC 67 12 has a strikingly unusual mass func- 
tion for stars below the turn-off ( Andreuzzi et alJl2001l) and 
this can only be expla ined if the cluster has lost more th an 
99% of its initial mass (Takah ashi & Portegies Zwari2 000) 4 . 

4. DISCUSSION 

For Galactic clusters, our estimates indicate that it is quite 
possible for all observed UCXBs to have been formed through 
RG-NS collisions (Table 1). For extragalactic clusters, we 
can crudely estimate the expected total number of UCXBs in 
a galaxy by integrating over the cluster mass function and 
assuming some average formation rate per NS in all clus- 
ters. As an example, consider the case of M87 (Jord an et alJ 
2004 ). We a dopt a powe r-law cluster IMF with a = -2 
( Kravtsov & Gnedin 2004). With an average formation rate 
K-ucxb ~ 2 x 10" 12 -4 x 10" n yr _1 per NS (assuming that 
M87 clusters have structural parameters distributed roughly 
like those of Galactic clusters) and fucxB ~ 10 6 yr (corre- 
sponding to L x > 10 37 ergs _1 , near the detectability limit for 
M87), we find that - 10- 100 UCXBs are expected in the 
1688 identified clusters, in rough agreement with the 58 de- 
tected LMXBs associated with these clusters. 

In several galaxies, the probability of finding a bright 
LMXB in a cluster appears to correlate strongly with clus - 
ter metallicity (iKundu et alJl2002l 12001 Uordin et al.ll2004l) . 
As there is no strong dependence of the RG-NS collision rate 
on metallicity, in our scenario, we have to interpret this trend 
as due to other factors, such as the metallicity dependence 
of the IMF, of the number of NSs formed in the cluster, or 
of the NS retention fraction. In other words, the metallicity 
dependence must appear through the number of NSs rather 
than through T^ucxb- This is difficult to verify, since there are 
no well established theoretical or observational predictions on 
how the IMF and NS natal kicks change with metallicity. Al- 
ternatively, a strong metallicity dependence of ?ucxb is also 
possible, with higher metallicity systems having longer life- 
times as bright sources (cf. Maccarone et al. 2004). 

Our conclusions are fairly robust, independent of assump- 
tions and in spite of some large theoretical uncertainties. We 
now examine a few of the most important ones. Binaries 
were neglected in our analysis. This implies that our esti- 
mated collision rate is a lower limit, as interactions involving 
binaries always increase this rate (Fregeau et al. 2004). How- 
ever, the effects of binaries on collision rates in very dense 
clusters today are likely to be small because core binary frac- 
tions in these clusters are very small, typica lly a few percent 
at most. This is known observationally ( Cool & Boltonl2002l) 
and expected theoretically (Fregeau et al. 2003; Ivanova et al. 
2004b). Another important assumption we made is that post- 
collision binaries do not circularize. As seen in Fig. 1, high 

4 After significant mass loss, and depending on its initial den- 
sity profile, the cluster could undergo strong gravothermal oscillations 
i Takahashi & Porteaies Zwarl 2000), so that a UCXB could also have formed 
when the core had a much higher density during a recent, brief episode of core 
collapse. 
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eccentricities are an important factor in keeping merger times 
short. However, one can also see directly from Fig. 1 that, 
even if all binaries were able to circularize quickly (compared 
to the GR merger time), a large fraction of post-collision sys- 
tems would still merge in less than the cluster age. Based on 
the results of Sec. 2 and the relation between post-collision 
semimajor axis and collision parameters derived from our 
SPH simulations, we estimate this fraction to be about 70%. 
Thus, even under the extreme assumption that all systems cir- 
cularize, the rate of UCXB formation would still be within 
a factor of 2 of the total RG-NS collision rate. One possi- 
ble further complication could come from the residual gas left 
bound to the RG core. All our collision calculations suggest 
that the mass left bound to the RG core is ~ 0.1M Q . Al- 
though there are many theoretical uncertainties, it is possible 
that this is sufficient to reconstitute a RG envelope (Castellani 
et al. 1994). In this case, the orbit would likely circularize, 



and stable MT from the reconstituted RG onto the NS would 
occur. However, the Roche lobe in the post-collision binary 
is smaller than the equilibrium radius of the RG, so that the 
MT proceeds on a thermal timescale and the corresponding 
bright LMXB phase lasts only ~ 10 5 yr, making detection un- 
likely. In addition, the total mass accreted by the NS will be 
only ~ 10~ 3 M Q , which is not sufficient to produce a recycled 
millisecond pulsar. 



ACKNOWLEDGMENTS 

We thank R. Bi, S. Fleming, V. Kalogera, M. Rosenfeld 
and B. Willems for helpful discussions. This work was 
supported by NSF Grants AST-0206276 and AST-0353997, 
NASA Grants NAG5- 12044 and NNG04G176G, and a Chan- 
dra Theory grant. 



REFERENCES 



Andreuzzi, G., De Marchi, G., Ferraro, F. R., Paresce, F., Pulone, L., & 

Buonanno, R. 2001, A&A, 372, 851 
Bildsten, L., & Deloye, C.J. 2004, ApJ, 607, LI 19 
Castellani, V., Luridiana, V., & Romaniello, M. 1994, ApJ, 428, 633 
Chakrabarty, D. 2004, to appear in Binary Radio Pulsars, ASP Conf. Ser., ed. 

FA. Rasio & I.H. Stairs | astro-ph/0408004 
Clark, G.W. 1975, ApJ, 199, CT41 

Cool, A.M., & Bolton, A.S. 2002, ASP Conf. Series, Vol. 263: Stellar 

Collisions, Mergers and their Consequences, 163 
Davies, M.B., Benz, W., & Hills, J.G. 1992, ApJ, 401, 246 
Deloye, C.J., & Bildsten, L. 2003 ApJ, 598, 1217 

Djorgovski, S. 1993, ASP Conf. Ser. 50: Structure and Dynamics of Globular 
Clusters, 373 

Deutsch, E.W., Margon, B., & Anderson, S.F. 2000, ApJ, 530, L21 
Fregeau, J.M., Cheung, P., Portegies Zwart, S.F, & Rasio, FA. 2004, 
MNRAS, 352, 1 

Fregeau, J.M., Giirkan, M.A., Joshi, K.J., & Rasio, FA. 2003, ApJ, 593, 772 
Gnedin, O.Y., Zhao, H., Pringle, J.E., Fall, S.M., Livio, M., & 

Meylan, G. 2002, ApJ, 568, L23 (supplementary online data at www- 

int.stsci.edu/ ognedin/gc/vesc.dat) 
Gnedin, O.Y., Hernquist, L., & Ostriker, J.P. 1999, ApJ, 514, 109 
Grebenev, S.A., et al. 2004, ATEL 350 

Hurley, J.R., Pols, O.R., & Tout, C.A. 2000, MNRAS, 315, 543 
Ivanova, N., Belczynski, K., Kalogera, V., Rasio, F.A., & Taam, R.E. 2003, 
ApJ, 592, 475 

Ivanova, N., Fregeau, J.M., & Rasio, F.A. 2004a, to appear in Binary Radio 
Pulsars, ASP Conf. Series, ed. F.A. Rasio & I.H. Stairs I astro-ph/0405382 i 



Ivanova, N., Belczynski, K, Fregeau, J.M., & Rasio, F.A. 2004b, MNRAS, 
in press 

Jordan, A., et al. 2004, ApJ, 613, 279 

Joshi, K.J., Nave, CP, & Rasio, F.A. 2001, ApJ, 550, 691 

Kalogera, V., & Webbink, R.F. 1996, ApJ, 458, 301 

Kalogera, V., Henninger, M., Ivanova, N., & King, A.R. 2004, ApJ, 603, L41 

Kroupa, P. 2002, Science, 295, 82 

Kravtsov, A.V., & Gnedin, O.Y. 2004, ApJ, in press 

Kundu, A., Maccarone, T. J., & Zepf, S. E. 2002, ApJ, 574, L5 

Kundu, A., Maccarone, T. J., Zepf, S. E., & Puzia, T. H. 2003, ApJ, 589, L81 

Lombardi, J.C., Dooley, K., Proulx, Z., Ivanova, N., Rasio, F. 2004, in prep. 

Maccarone, T.J., Kundu, A., & Zepf, S.E. 2004, ApJ, 606, 430 

Peters, PC. 1964, Phys. Rev., 136, 1224 

Podsiadlowski, P., Rappaport, S., & Pfahl, E.D. 2002, ApJ, 565, 1107 
Pryor, C. & Meylan, G. 1993, ASP Conf. Ser. 50: Structure and Dynamics 

of Globular Clusters, 357 
Pylyser, E., & Savonije, G.J. 1988, A&A, 191, 57 
Rasio, FA., Pfahl, E.D., & Rappaport, S. 2000, ApJ, 532, L47 
Rasio, FA., & Shapiro, S.L. 1991, ApJ, 377, 559 
Richer, H.B., et al. 2004, AJ, 127, 2771 

Richtler, T, Grebel, E.K., & Seggewiss, W. 1994, A&A, 290, 412 
Takahashi, K. & Portegies Zwart, S. F. 2000, ApJ, 535, 759 
van der Sluys, M. V., Verbunt, F, & Pols, O. R. 2005, A&A, 431, 647 
Verbunt, F. 1987, ApJ, 312,L23 

Webbink, R. F. 1985, IAU Symp. 113: Dynamics of Star Clusters, 113, 541 



Note added in proof. - Bildsten & Deloye (2004) pointed out that the observed break around 5 x 10 38 ergs _I in the X-ray 
luminosity functions of elliptical galaxies could be explained naturally if all donors in UCXBs had initial masses clustered near 
~ 0.1 Mq. This is precisely what is predicted by our scenario. 



